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The fluctuation conductivity <r s in bulk superconductors with non s-wave pairing and with 
nonmagnetic disorder of strength D is studied at low T and within the Gaussian approximation. 
It is shown by assuming a quasi two-dimensional (2D) electronic state that, only if the gap 
function d^(p) is, as in a 2D p-wave pairing state, linear in the in-plane (relative) momentum 
p^, the in-plane fluctuation conductivity on the line D = D c is weakly divergent in low T limit. 
The present result may be useful in clarifying the true gap function of spin-triplet Si'2Ru04 
through resistivity measurements. 
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Conventionally, the pairing state of a superconducting (SC) phase is identified by examining 
microscopic properties deep in the SC phase or by directly finding a pairing mechanism consistent 
with the electronic state of the material. In contrast, it is usually difficult to obtain information on 
the pairing state through physical properties near a SC transition such as the fluctuation effects. For 
example, the fluctuation conductivity cr s accompanying a thermal transition depends only on the 
dimensionality© of fluctuation except for a constant numerical prefactor of its expression, which is 
noti'& easy to examine experimentally with high precision. But how about near a SC transition at 
zero temperature (T = 0) ? From this point of view, we examine the nature of o~ s in unconventional 
(non s-wave) superconductors near a quantum critical point, (D,T) = (D c ,0), which is the end 
point of disorder-driven transition line T C (D), where D is the strength of microscopic (nonmagnetic) 
disorder. Particular attention is paid to a s on the specific line D = D c in the D — T phase diagram. 
We find in the Gaussian approximation that, only in the case of pairing of p-wave typejl' <r s on the 
D = D c line is weakly divergent in low T limit. 

This work was motivated by the unresolved issue of the pairing symmetry of Sr2RuO4.II) In 
spite of various proposals, the pairing state of this material has not been determined thus far. 
A key issue is to determine where on the Fermi surface the nodes of SC energy gap are present. 
Having Sr2Ru04 in mind, we assume that the electronic states are of quasi 2D character, while the 
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dimensionality of low energy modes such as the superconducting fluctuation is three-dimensional 
(3D). Further, all calculations will be carried out within the weak coupling approximation, and the 
effects of disorder on the fluctuation, as well as on the normal properties, will be treated within 
the Born approximation since the nature of fluctuation is determined by the electronic details in 
the normal phase. A key fact in the present analysis is that, when the gap function is of p-wave 
type (according to the classification defined in ref. 4), a part of the lowest order gradient terms 
in the Gaussian Ginzburg — Landau (GL) action is proportional to a ln(l/Tr)-factor arising from 
the diffusion propagator induced by the disorder, where r is the lifetime of quasiparticles, while for 
gap symmetries with higher harmonics (i.e., cf-wave, /-wave types^ or others) the corresponding 
ln(l/Tr)-factor appears only in higher gradient terms that are unnecessary for obtaining the leading 
term of a s . 

As a microscopic model, we start with the BCS Hamiltonian of a layered material with a random 
potential term 

T J) I J 

-\Y. E Wp, pH(p + WA-p-W'(-p'-W(p' + 1 (i) 

<y,<y' p,p'q 

where a a (r±, I) is the annihilation operator of electron with spin a and at the in-plane coordinate 
r± on the l-th plane, the attractive interaction Vbcs(p> p') is assumed to be, as usual, separable like 
\g\w(p)w(p'), and p ± implies p±q/2. The random potential ui(r±) obeys the Gaussian ensemble: 
ui(i'±) = 0;ui(r±)ui'(r'j_) = (2irN2(0)T)~ 1 5(r± — r' ± )5i : i', where the overbar denotes the random 
average, and No(0) is the 2D density of states at the Fermi surface. The impurity scattering was 
assnmeci at this stage to b e isotropic. We nave verified that an inclusion of p-wave scattering ratj 
t^ 1 does not affect the main features in eq. @ becoming relevant to our derivation of a s (D = D c ). 
The random-averaged Green's function is, as usual, given by 

G(p,£n) = -. n , , n i 1 7-, (2) 

ie n (l + (2r|e n |) x ) - i p 
with a fermionic Matsubara frequency e n . The dispersion in eq.(p|) takes the quasi 2D form 

£, P = £ P± - E F + tcos(p z s), (3) 

where s is the inter layer distance and p = (p±,p z )- Throughout this paper we take units Ti = = 
1. Also we neglect the localization effect on noninteracting electrons and treat the disorder effect 
within the Born approximation. 

To illustrate the points essential to our results, let us first explain the analysis of a case of d-wave 
type defined by Vbcs = \g\w d {p)w* d {p') with w d (p) = (p 2 x -p 2 y )h d (p z s). Here, h d (p z s) is assumed to 
satisfy {h d ) = s f h d {p z s) > so that h d may be p^-independent, reflecting the quasi 2D dispersion 
©• 
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Fig. 1. Feynman diagrams defining the renormalized vertex Wq(p) (the hatched region), where the double line, 
the solid line with arrow, and the dashed line denote, respectively, the pair-field, the electron propagator, and the 
impurity line carrying (2tttN2(0))~ 1 . The first term of r.h.s. implies the bare vertex Wd(p)- 



Following the standard procedure^ to decoupe the BCS interaction in terms of the pair-field 
which is related to the ci-vector as d^(p) = u>d(p)A M , where fi is the spin index, the Gaussian GL 
action takes the form 

* S = EE /C_1 (^^)|A M (Q,^)| 2 , (4) 
where K~ l = flglA^O)) -1 - ir(q,u v ), and 

ir(q,u v ) =T£ J de p± j ^l w *(p) Wq (p)g(p + ,s n+ )g(-p_,-s n ^, (5) 

where the renormalized vertex Wq(p) is defined in terms of the bare vertex Wd(p) as indicated 
byS* Fig. 1. In general, features of crystal structure appear in the dispersion e p± which may make 
the derivations of eq.(|5|) and Wq(p) complicated. If an isotropic Fermi surface is assumed, for 
simplicity, Wq(p) is given by 



2 /2 

W q (p) = w d {p) - ^-{h d ){qt - q 2 y ) T q (2e n + u v ), (6) 



where Tq{2e n + uj v ) = (r|2e n + uj u \ + l 2 q\/2 + 2t 2 r 2 sin 2 (g 2 s/2)) _1 is the diffusion propagator, and 
uj u is a bosonic Matsubara frequency. Since the q-dependence in the second term [i.e., q 2 . — g 2 in 
eq. (||) ] is essentially determined by the p-dependence of w d , the effect of crystal structure is, 
except for its role leading to the form of Wd(p), unimportant in deriving the renormalized part 
corresponding to the second term of eq.(||). An effect of crystal stucture on the GL gradient terms 
unaccompanied by the diffusion propagator will be commented on later together with an effect of 
the GL quartic term neglected here. 

It is easily seen that the last term of eq. leads to higher gradient terms like ~ (g 2 — g 2 ) 2 in 
GL action. A (^-dependence neglected above, which may appear in (h^) of Wq, is reflected only in 
a still higher gradient term like (g 2 — q^) 2 q\. Since such higher gradient terms are irrelevant to the 
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leading terms of interest in fluctuation conductivities, we can omit the last term of eq.(^) at this 
stage. Then, JC(q, uj u ) simply becomes 

KT\q,u} v ) = /i + 7iM + fo?± + 4 (^f) sm 2 {q z s/2). (7) 

Here /x = 2tt 2 (Tt) 2 /3 + (D - D c ) /D c with D/D c = 1.78/(ttT c0 t), 71 ~ r, $ = l 2 /2 with 
I = vpr, £ 2 = (trs) 2 /2, and T c o is T c in clean limit (D = 0). Note that the fact J p Wd(p) = is 
essential for obtaining the T c -reduction induced by disorder. Hereafter, we only take the dissipative 
frequency term& into account. Inclusion of other dynamical terms is found not to change our main 
conclusions. The above-mentioned feature that the contribution accompanied by the diffusion 
propagator is reflected only in higher gradient terms of GL action is common to all pairing states 
with harmonics higher (angular momentum) than those of p-wave type,& and hence, the results 
from eq. (9), below, for the case of cf-wave type are also valid for other states with higher harmonics 
including states of /-wave type. 

Next, let us calculate the fluctuation conductivity cj s j_ for an in-plane current using the Gaussian 
GL action and the Kubo formula. As is well known, this is nothing but the Aslamasov — Larkin 
(AL) fluctuation conductivity & After the analytical continuation, it is expressed by 

2tt$ f 1 f°° d 2 q ± 2 r/ d dq z r R i2 

16 ^ T tfM { w))i** - , i, ( 8) 



" 3^00 ' ' "J (2VT) 3 [(l + g 2)2 +w 2]2' 

where Rq = 7rh/2e 2 = 6.45 (kf2), K R (q,u>) = [/C _1 (g, 0) — i7io;] _1 is the retarded fluctuation 
propagator, and, in moving to the last line, a 3D-type behavior of fluctuation was assumed. This 
assumption is valid when (0 <) y / /Io < 2£ c o/s. The weight of quantum fluctuation in cr s -expressions 
is measured by the function f(x) = x 2 / sinh 2 (a;). 

First, let us comment on the strong disorder side D > D c in which ^0 remains positive in T — > 
limit. It is easily seen that o" s j_ vanishes like© T 2 in T — > limit in this case. This is the simplest 
example of the insulating^ fluctuation conductivity at T = 0. On the other hand, on the weak 
disorder side in which D < D c , fiQ vanishes at a positive temperature T C (D) on cooling. Near T C (D) 
where /io *C 2 r yiT c (D), the familiar divergent (superconducting) behavior!) a s ±_ ~ (T — T C (D)) -1 / 2 
is found. Now, let us consider a s ± on the line D = D c where fiQ = 2tt 2 (Tt) 2 /3. The ratio fio/2jiT 
in this case also tends to vanish on cooling and becomes less than unity in T < T C r — 0.2T c o, 
and thus the same classical limit (f(x) — > 1) as in D < D c dominates at D = D c in low T limit. 
Through a detailed analysis, we find the remaining quantum correction in <t s j_ to be smaller by 



the factor ~ ^T/T c q. Further, the condition of the 3D approximation, yfilo < 2£ c o/s, is rewritten 
in this case as T < — s/3t/{2ir). Then, in T <C min(T C r 1 ' > , cr s ±(D = D c ) approaches a 
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T-independent value 



a s± (D = D c ) 



V3 T ( 
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4vrV 2£ c0 (D = D c ) 7R Q s t 
on cooling. By combining this with the D ^ D c results mentioned above, one will notice that the 
situation around D = D c in the present 3D case is similar to that of the 2D insulator-superconductor 
transition® A similar result is obtained for the out-of-plane conductivity cr s || of which the value 
at D = D c is given by 



a s]] (D = D c ) =a, 



* £c0 
6) 
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a s± (D = D c 



(10) 



Now, we consider the cases of pairing of p-wave type.tr Let us assume the attractive interaction to 
have a general form V B cs(p,p') = \g\ Y l j=x, y w j(p) w j(p')> where vij = Pjhj(p z s) with J Pz hj(p z s) > 
0. Then, the d-vector has the form d^ip) = J2j=x,yPj^j(Pz s )^^,j- This includes the candidates 
for Sr2Ru04 (triplet) pairing state,0'@) d^{p) = J2 P =±i A p (p x + \pp y ){c — cosp z s). As in this 
example, we will assume \hj(p z s)\ to be independent of j for convenience of presentation. Under 
these apparatus, it is elementary to obtain the corresponding expressions to eqs. (||) and @ again 
in terms of the isotropic dispersion. The resulting Gaussian action is merely given here: 



S (p) ^ ^2(0) 



where 



K T x {q,u3 v ) 



£c0 



Kl\q, u v ) = K T \q, u v ) + $C(q 2 ± ,T) q\, (12) 

cjj = qj/\q±\, and C{q\,T) = 2(1 + ln(47rrT + q^l 2 )^ 1 ). The enhancement factor C{q\,T) arising 
from the diffusion propagator appears only in the longitudinal part of the lowest order terms in the 
gradient because the term corresponding to the second term of eq. @ is vectorial and proportional 
to q. An additional gradient term which may appear depending on the form of e p± will be discussed 
later. The in-plane conductivity a s ± in this case becomes 



<r$(P = D c ) =3x 



27t£ 4 f d 3 q 



q\ 



RqTJ (2vr) 3 J u sinh 2 ($,) 
(C{ql,T)f 



(lmK^(q,u;)) 2 + (1 + C(q 2 ± ,T)) 2 (lmK?(q, 



(lmK*(q,Lu)) (hnKj^q, 



D=D C 



6cj* (lnC(T) + const.) 



(13) 
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on the line D = D c and in low T limit, where C(T) = C(q\ = 0, T) ~ ln(rT) -1 , and the factor 3 
in the first line is due to spin degeneracy in the spin-triplet channel. In eq. (13), we have neglected 
contributions associated with an additional term of supercurrent arising from the g_|_-dependence 
of C{q\,T), which merely become, at most, of the order ct* C -1 (T) In C(T). On the other hand, 
the out-of-plane conductivity is given by 



a^(D = D c )^3a:(l + C' 1 (T))^ 



(14) 

D=D C 



(p) 

The low T behaviors of ak FJ, s in D 7^ D c are similar to those in the case of non p-wave type: 
ai p) (D > D c ) -» T 2 , while o£°(D < D c ) ~ (T - T c (D))~ l l 2 when T -» T C (L>)+. As eq. (13) shows, 
however, on the line D = D c is, in contrast to eq. (9), weakly divergent (superconducting) in 
T —> limit. 

We should note that it is the sum of a fluctuation contribution cr s and a normal one <r n which is 
directly measured in experiments of the type performed in ref. (13). If the electronic structure is of 
3D-type, the residual normal conductivity is of the order R^Eprkp, implying that the ratio cr*/a n 
will be small, 0((T c o/-Ef) 2 )> f° r 3D-type Fermi surfaces. In the present case with a quasi 2D Fermi 
surface, however, this ratio is larger: since a n ~ Rg 1 EpT/2s, we have c7*/<7 n ~ y/Es / (8tt Eft £co(D = 
D c )), which becomes 0.67rT c 2 /(£ , Fi) in the present model. In a system with £ c0 ~ s (i.e., t ~ T c0 ) 
such as Sr2Ru04, this ratio is 0(T c o/-Ef)- Although the ratio T c o/Ep in Sr2Ru04 is small (~ 10~ 2 ), 
it will be possible to separate a B from <r n by applying a magnetic field as long as a n is insensitive to 
the field, because the fluctuation contribution er s , in particular near D = D c , will be easily erased 
by applying a low magnetic field. 

A few comments on the above derivation of <t s are in order. First, two main effects of crystalline 
anisotropy were neglected above. As one of them, the p-dependence of the gap function dn(p) was 
linearized above, and, as an example, Wj(p) should be proportional not to pj but to s'm pj in a 
real system. However, it is clear that the gradient terms, in particular its quadratic terms, of GL 
action do not essentially change due to this replacement. Next, the other gradient term which may 
appear in eq. (11) as a result of reduction of space symmetry should also be examined, which is 
typically given by 

\8 X A^ X \ 2 + IdyA^J 2 (15) 

in zero field. The resulting GL gradient terms have three independent coefficients and hence, 
are essentially the same as those in the so-called 2D representation© as long as the gauge field is 
spatially constant (note that a conductance for a uniform current in zero magnetic field is defined in 
terms of a spatially constant gauge field). However, the inclusion of eq. (|i"5|) does not qualitatively 
change our main results (9) and (13). A key point for reaching this conclusion is to notice that 
the contribution accompanied by the diffusion propagator leads not to the term (15) but only to 
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a longitudinal term ~ |divA| 2 and that, as in clean limit, the term (15) can arise simply as a 
combined result of crystal anisotropy in the dispersion e p± and of the Wj(p)-form. Then, it is 
straightforward to see, by power counting in q-integrals, the appearance in cases of p-wave type of 
a lnln(Tr) _1 -divergence like that in eq. (13) even when the term fli~5| ) is present. 

We have focused above on the AL term of Gaussian fluctuation conductivity. It is not difficult to 
examine other contributions, Maki — Thompson^ term <5<tmt an d the so-called DOS term 5o"dos- 
We find that, in T — ► limit, they vanish on the D = D c line in the manners 5<7mt ~ Tt and 
^cdos ~ (^ T ) 2 ) an d thus, we can neglect them as far as the behavior in low T limit is concerned. 

Fluctuation conductivity near the disorder-induced quantum critical point was also examined in 
ref.15, where the p-wave case does not seem to have been considered. The corresponding result to 
eq. (9) was concluded in ref. 15 to be rather proportional to T 1 / 4 . This difference originates from 
the assumption implicit in ref. 15 that the mass renormalization 5fi ~ T D I 2 in D-dimensions in 
the one-loop order0) is larger than the bare one fio ~ T 2 . However, the former arises from the 
quartic (nonGaussian) term of GL action. In fact, it is verified that the above assumption 5jjL > fiQ 
is vaiid in T < r c0 eg>, where «g> a is the 3D Ginzburg number ™i„g the 

temperature width AT/T c q of the thermal critical region. Namely, if this temperature width of 
O(T c 3 /i?p) is negligible as in Sr2Ru04, the T^-behavior of a s (D = D c ) will not be observable at 
accessible temperatures, but rather, the present Gaussian results should be observed at low T. 

We emphasize that the result, eq. (13), in cases of p-wave type is valid only when the Fermi 
surface and thus, the gap function d(p) are of quasi 2D-type. If the Fermi surface and hence, 
d(p) is essentially 3D-type (e.g., d(p) = J2j=x,y,z ^jPj)i a s{D = D c ; T — » 0) is nondivergent 
and is not useful in distinguishing the pairing states. By contrast, if d(p) is ID-type and, as its 
example, is approximated by sin(p a .), the fluctuation conductivity cr sx (D = D c ) in the x-direction is 
divergent like (h^l/Tr)) 1 / 2 , while the corresponding ones in the y and z directions instead vanish 
like (Intl/Tr))- 1 / 2 in low T limit. 

In conclusion, the fluctuation conductivities near a disorder-induced quantum critical point of 
a non s-wave bulk superconductor with quasi 2D-like electronic structure were examined in the 
Gaussian approximation and in zero field. When the dependence of the gap function on the in- 
plane momentum p^ is of a higher order, as in cases of c?-wave or /-wave symmetry, a s at the 
crtitical value of disorder approaches a finite value in low T limit, just like the 2D nonGaussian 
result ,0 while it is weakly divergent on cooling when the gap function is linear in p ± or sin(p^) 
(of a p-wave type). The present result may be useful for clarifying the form of pairing function of a 
material with spin-triplet pairing and a low dimensional Fermi surface and, for example, clarifying 
which of the /-wave gap function0 oc p x db ip y and the function© oc p x Py(Px ± Wy) ls favorable as 
the pairing state of Sr2Ru04. 
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